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Abstract —  In  this  paper  we  introduce  the  notion  of  optimiza¬ 
tion  under  control  and  communication  constraint  in  a  robotic 
network.  Starting  from  a  general  setup,  we  focus  our  attention 
on  the  problem  of  achieving  rendezvous  in  minimum  time 
for  a  network  of  first  order  agents  with  bounded  inputs  and 
limited  range  communication.  We  propose  two  dynamic  control 
and  communication  laws.  These  laws  are  based  on  consensus 
algorithms  for  distributed  computation  of  the  minimal  enclosing 
ball  and  orthotope  of  a  set  of  points.  We  prove  that  these  control 
laws  converge  to  the  optimal  solution  of  the  centralized  problem 
(i.e.,  when  no  communication  constrains  are  enforced)  as  the 
bound  on  the  control  input  goes  to  zero.  Moreover,  we  give  a 
bound  for  the  time  complexity  of  one  of  the  two  laws. 

I.  Introduction 

The  interesting  aspect  of  motion  coordination  consists  in 
combining  together  problems  from  control  and  communica¬ 
tion  theory.  The  main  difficulty  deals  with  integrating  the 
sensing,  computing,  communication  and  control  aspects  of 
problems  involving  groups  of  mobile  agents.  A  well  known 
problem  in  control  theory  is  optimal  control.  Roughly  speak¬ 
ing,  it  consists  in  finding  a  feedback  law  that  minimizes  some 
cost  functional  under  some  inputs  and  dynamics  constraint. 
In  this  paper  we  introduce  the  notion  of  optimal  control  and 
communication  for  a  network  of  robotic  agents.  We  want 
to  study  how  to  solve  an  optimization  problem,  in  presence 
of  both  the  usual  motion  constraints  and  the  communication 
ones.  In  particular  this  paper  is  a  preliminary  contribution 
towards  what  might  be  loosely  referred  to  as  “distributed 
geometric  optimization.”  In  fact  many  optimization  problems 
for  robotic  networks  can  be  shown  to  be  equivalent  to  the 
computation  of  geometric  shapes.  While  in  a  centralized 
setting  the  solution  is  usually  simple,  the  problem  becomes 
very  complicated  when  it  must  be  solved  in  a  distributed 
way.  Distributed  computation  over  network  has  been  largely 
studied  for  fixed  topologies;  e.g.,  see  [1], 

In  this  paper  we  point  our  attention  on  the  well  known 
rendezvous  coordination  task  and  look  for  solutions  that 
solve  such  task  in  minimum  time.  We  look  for  distributed 
solutions  in  networks  of  mobile  agents  with  first  order 
dynamics,  bounded  inputs  and  limited-range  communication. 

The  “multi-agent  rendezvous”  problem  and  a  “circumcen- 
ter  algorithm”  have  been  introduced  by  Ando  and  coworkers 
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in  [2].  The  algorithm  proposed  in  [2]  has  been  extended  to 
various  synchronous  and  asynchronous  stop-and-go  strate¬ 
gies  in  [3].  A  related  algorithm,  in  which  connectivity  con¬ 
straints  are  not  imposed,  is  proposed  in  [4].  In  [5]  the  class 
of  “circumcenter  algorithms”  has  been  studied  in  networks 
of  agents  whose  state  space  is  Rd,  for  arbitrary  d,  and  with 
communication  topology  characterized  by  proximity  graphs 
spatially  distributed  over  the  disk  graph.  In  [6]  the  (time  and 
communication)  complexity  of  this  and  other  algorithms  has 
been  studied.  All  these  coordination  schemes  are  memoryless 
(static  feedback).  In  this  paper  we  want  to  explore  dynamic 
control  and  communication  laws  in  order  to  approximate 
the  optimal  solution  of  the  minimum  time  rendezvous.  In 
particular  the  control  and  communication  laws  is  based 
on  reaching  consensus  on  some  logic  variables  and  at  the 
same  time  moving  toward  the  current  estimation.  A  similar 
approach  was  used  in  [7]  where  the  agents  try  to  reach  a 
consensus  on  a  set  of  variables  called  coordination  variables. 

Studying  the  minimum  time  rendezvous  problem  in  the 
centralized  setting  we  show  that,  depending  on  the  norm 
used  to  bound  the  control  input,  the  optimal  solution  consists 
of  moving  toward  the  center  of  the  minimal  enclosing  ball 
(bound  on  L2  norm)  and  toward  the  center  of  the  minimal 
enclosing  orthotope  (bound  on  the  infinity  norm)  of  the 
points  located  at  the  initial  position  of  the  agents. 

Our  main  result  is  the  design  of  a  control  and  com¬ 
munication  law  based  on  a  consensus  algorithm  for  the 
distributed  computation  of  the  minimal  enclosing  ball  and 
the  minimal  enclosing  ortho  tope  of  a  set  of  points.  We  prove 
the  correctness  of  the  two  consensus  algorithms  and  provide 
a  bound  on  the  time  of  convergence  for  the  orthotope  case. 
Then  we  prove  that  the  control  and  communication  law  that 
combines  the  consensus  with  the  motion  law  converges  to  the 
optimal  solution  as  the  control  bound  goes  to  zero.  Moreover, 
for  the  problem  with  input  bounded  by  the  infinity  norm 
(corresponding  to  the  computation  of  the  minimal  enclosing 
orthotope),  we  prove  that  the  control  and  communication  law 
is  a  constant  factor  approximation  of  the  centralized  optimal 
solution. 

In  Section  m  we  introduce  a  formal  model  of  robotic 
network  inspired  by  the  one  introduced  in  [6].  Moreover, 
we  define  the  optimal  control  and  communication  problem. 
In  Section  [TO]  we  characterize  the  solution  of  minimum  time 
rendezvous  in  a  centralized  setting.  In  Section  HV1  we  define 
the  FloodMEB  and  FloodMEO  algorithms  for  the  distributed 
computation  of  minimal  enclosing  ball  and  orthotope,  prove 
their  correctness  and  give  bounds  on  time  complexity.  Sec¬ 
tion  ^contains  the  control  and  communication  laws  based  on 
the  consensus  algorithms  described  in  Section  EY]  Finally  in 


Section  IV-CI  and  Section  ED  we  show  simulations  and  draw 
the  conclusions  with  future  perspectives. 

II.  Preliminary  developments 

In  this  section  we  recall  the  concepts  of  network  of 
robotic  agents,  coordination  tasks  and  complexity  measures, 
and  introduce  the  notion  of  optimization  under  motion  and 
communication  constraints. 

A.  Notation 

We  let  N,  No,  and  K+  denote  the  natural  numbers,  the 
non-negative  integer  numbers,  and  the  positive  real  numbers, 
respectively.  We  let  JI  r e { 1  n}  denote  the  Cartesian 
product  of  sets  Si, ... ,  Sn.  For  pel,  we  let  |_ p\  and  |"p~| 
denote  the  floor  and  ceil  of  p.  For  r  £  K+  and  p  £  Rd, 
we  let  B(p,r)  denote  the  closed  ball  centered  at  p  with 
radius  r,  i.e.,  B(p,r)  =  {q  £  ffi.d  |  \\p  —  q H2  <  r}  and 
C(p,r )  denote  the  closed  hypercube  centered  at  p  with 
sides  of  length  r  and  parallel  to  the  coordinate  axes,  i.e., 
C(p,  r)  =  {q£Rd  |  ||p  -  q||oo  <  r}. 

For  /,  g  :  N  — >  R,  we  say  that  /  £  0(g)  (respectively, 
/  £  Cl(g))  if  there  exist  uq  £  N  and  k  £  R+  such  that 
|/(n)|  <  k\g(n)\  for  all  n  >  no  (respectively,  |/(n)|  > 
k\g(n)\  for  all  n  >  no).  If  f  £  O(g)  and  /  £  f 1(g),  then  we 
use  the  notation  /  £  0(p). 

Next,  we  briefly  review  some  useful  proximity  graphs. 
Given  rcmm  (E  M+,  the  disk  graph  (/disk(rCmm),  respec¬ 
tively  cube  graph  ^cube(fcmm),  is  the  state  dependent 
graph  on  W1  defined  by  the  following  statement:  for  any 
pointset  {pM, . . .  ,p^}  C  Rd,  the  pair  ( i,j )  is  an  edge 
in  f/diskOcmm)  •  ({p[1], . . .  ,p[n]}),  respectively  f7cube(rCmm)  ■ 
({pi1!, . . .  ,p!nl}),  if  and  only  if  i  7^  j  and 

||pM  -  P[j]  || 2  <  rcmm  pw  -  p[j]  £  B(0d,  rcmm), 

respectively 

||PW  -V[3]  II 00  <  Amm  PW  -P[jl  £  C(Od,rcmm). 

Another  useful  graph  is  the  complete  graph  (?cmpi,  i.e.,  the 
graph  with  edges  between  any  pair  of  nodes. 

Finally,  given  a  graph  Q  (even  not  state  dependent),  we 
denote  with  distg(j,j)  the  topological  distance  between  i 
and  j,  i.e.,  the  minimum  number  of  agents  to  go  from  i  to  j 
in  the  graph  Q.  We  define  diamg,  the  diameter  of  Q,  to  be 
the  maximum  topological  distance,  dist q(i,j),  for  all  (i,  j). 

B.  Modeling  a  network  of  robotic  agents 

We  describe  a  (uniform)  network  of  robotic  agents  using 
the  formal  model  introduced  in  [6]  modified  for  the  discrete 
time  case.  The  network  is  modeled  as  a  tuple  (/,  A,  Bcmm). 
I  =  {l,...,n}  is  the  set  of  unique  identifiers  (UIDs); 
A  =  {A^}ieI  =  {(X,  U,X0,f)}iei  is  called  the  set  of 
physical  agents  and  is  a  set  of  control  systems  consisting  of 
a  differentiable  manifold  X  (state  space),  a  compact  subset 
U  of  Rm  (input  space),  a  subset  Xq  of  X  (set  of  allowable 
initial  states)  and  a  (sufficiently  smooth)  map  f  :  XxU  — >  X 
describing  the  dynamics  of  ith  agent;  A’cnn,  :  Xn  —*  I  x  I 
is  called  the  communication  edge  map. 

The  robotic  network  evolves  according  to  a  discrete-time 
communication  and  motion  model. 


Definition  2.1  (Control  and  communication  law):  Let  S 
be  a  robotic  network.  A  (uniform,  synchronous,  dynamic) 
control  and  communication  law  CC  for  S  consists  of  the 
sets: 

(i)  L,  a  set  containing  the  null  element,  called  the 
communication  language',  elements  of  L  are  called 
messages; 

(ii)  W,  set  of  values  of  some  logic  variables  i  £  I; 

(iii)  Wo  C  W,  subsets  of  allowable  initial  values; 
and  of  the  maps: 

(i)  msg  :  X  x  W  X  /  — »  L,  message-generation  function; 

(ii)  stf  :  W  x  Ln  — >  W,  called  state-transition  function; 

(iii)  ctl  :  X  x  W  x  Ln  — >  U,  called  control  function.  □ 

Roughly  speaking  this  definition  has  the  following  mean¬ 
ing:  for  all  i  £  /,  to  the  /th  physical  agent  corresponds  a 
logic  process,  labeled  i,  that  performs  the  following  actions. 
First,  at  each  communication  round  the  /th  logic  process 
sends  to  each  of  its  neighbors  in  the  communication  graph  a 
message  (possibly  the  null  message)  computed  by  applying 
the  message-generation  function  to  the  current  values  of 
and  w®.  After  a  negligible  period  of  time,  the  ith  logic 
process  resets  the  value  of  its  logic  variables  w  ^  by  applying 
the  state-transition  function  to  the  current  value  of  w®,  and 
to  the  messages  received  at  time  t.  Between  communication 
instants,  the  motion  of  the  ?th  agent  is  determined  by 
applying  the  control  function  to  the  current  value  of  x^l\  and 
the  current  value  of  u/M.  This  idea  is  formalized  as  follows. 

Definition  2.2  (Evolution  of  a  robotic  network ):  Let  S  be 
a  robotic  network  and  CC  be  a  control  and  communication 
law  for  S.  The  evolution  of  (S.  CC )  from  initial  conditions 

M  r 

xf  £  Xq  and  w0  £  Wo,  i  £  I,  is  the  set  of  curves  xm  : 
N  — ►  X  and  w/M  :  N  — >  W,  i  £  I,  satisfying 

(t  +  1)  =  /(*W (t),  ctl(*N (t ), (t  +  1),  (t))) , 

where,  for  i  £  I, 

wW(t  +  l)  =  stf{wW(t),yW(t)), 

with  the  conventions  that  cc^!  (to)  =  ^  and  tt)!*!(fo)  = 
i  £  I.  Here,  the  function  yW  :  N  — >  Ln  (describing  the 
messages  received  by  agent  i)  has  components 

msg(:rw(f),u/w(f),*),  if  (i,  j)  £  Ecmm,  Q 
null,  otherwise. 

In  the  paper  we  consider  the  following  network.  Each 
agent  i  occupies  a  location  pM  £  Rd,  d  £  N,  and  moves 
according  to  the  first  order  discrete-time  integrator 

pN(i  +  l)=pw(t)+u[i](t).  (1) 

The  communication  edge  map  can  be  either  the  one  arising 
according  to  the  disk  graph,  F^isk,  or  the  one  according 
to  the  cube  graph,  Ecu be.  Each  control  takes  values  in 
a  bounded  subset  of  Rrf,  that  can  be  either  f3(0,rctr)  or 
C(0,rctr),  i.e.,  ||zt!ll||2  <  rctr  or  ||mw||oo  <  rctr.  Notice  that, 
in  general,  the  type  of  communication  edge  map  and  the 
type  of  control  bound  are  not  related.  Finally  the  control 
and  communication  law  will  be  defined  depending  on  the 
coordination  task. 


C.  Coordination  tasks  and  time  complexity 

We  are  ready  to  define  the  notion  of  task  and  of  task 
achievement  by  a  robotic  network. 

Definition  2.3  ( Coordination  task):  Let  S  be  a  robotic 
network.  A  (static)  coordination  task  for  S  is  a  map  T  : 
Xn  — >  {true,  false}.  Additionally,  let  CC  a  control  and 
communication  law  for  S.  The  law  CC  achieves  the  task  Tif, 

M  M 

for  all  initial  conditions  x0  G  X0  and  w0  G  Wq,  iG/,  the 
corresponding  network  evolution  t  >— >  (x(t),w(t))  has  the 
property  that  there  exists  T  G  N  such  that  T\x(t))  =  true 
for  all  t>T.  □ 

We  are  finally  ready  to  define  the  notion  of  time  com¬ 
plexity  as  the  minimum  number  of  communication  rounds 
needed  by  the  agents  to  achieve  the  task  T with  CC. 

Definition  2.4  (Time  complexity):  Let  S  be  a  robotic  net¬ 
work  and  let  The  a  coordination  task  for  S.  Let  CC  be  a 
control  and  communication  law  for  S  compatible  with  T. 
The  time  complexity  to  achieve  T  with  CC  from  Xq  G  Xq  is 

TC(T, CC,xo)  =  inf  {T  G  N  |  T(x(t))  =  true  ,  Vi  >  T} 

where  1 1— >  (x(t),w(t))  is  the  evolution  of  (S,CC)  from  the 
initial  condition  (atojitfo)- 

The  time  complexity  to  achieve  T with  CC,  TC(7 ,CC),  is 
the  maximum  TC(T,  CC,  xq)  over  all  initial  conditions  xq. 

D.  Optimal  control  and  communication  in  robotic  networks 

Having  defined  a  coordination  task  for  a  robotic  network, 
we  can  ask  whether  such  task  can  be  accomplished  minimiz¬ 
ing  some  cost  functional.  In  what  follows  we  will  introduce 
the  notion  of  optimal  control  and  communication  problem 
and  of  optimal  control  and  communication  law  as  solution 
of  the  problem. 

Definition  2.5  (Optimal  control  and  communication): 
Given  a  task  T  and  a  cost  functional  J(u(-),x(T),T),  an 
optimal  control  and  communication  problem  is  the  following: 

minimize„(.)iX(0)lX(T),T  J(u(-),x{T),T) 

J(u(-),x(T),T)  =  Et=o(1(x(t)>u(t))  +Q(x(t))’ 

subj.  to 

(i)  (x(-),u(-))  is  an  input-state  trajectory  of  A, 

A  =  {A^}ieI\ 

(ii)  i  and  j  can  communicate  if  and  only  if 

(i,j)  G  Ecmm(x[1](t),...,x^(t)); 

(iii)  T(x(t))  =  true  for  all  t  >  T,  T  G  N. 

where  l  :  Xn  x  Un  — >  R.  is  a  sufficiently  smooth  and 
nonnegative-valued  function,  called  stage  cost,  and  g  : 
Xn  — >  M  has  the  same  properties  plus  g(x)  =  0  for  all 
x  G  Xn  such  that  D(x)  =  true  (for  an  admissible  CC).  □ 
We  say  that  a  control  and  communication  law  CC  is 
optimal  with  respect  to  the  coordination  task  T  and  the 
cost  functional  J,  if  it  solves  the  above  optimal  control  and 
communication  problem. 

We  call  CC  a  centralized  optimal  control  and  communica¬ 
tion  law  if  it  solves  the  optimization  problem  for  a  network  of 
robotic  agents  that  communicate  according  to  the  complete 
graph,  i.e.,  the  communication  edge  map  is  Ecmp\. 


Remark  2.6:  The  centralized  solution  of  an  optimal  con¬ 
trol  and  communication  problem  is  the  classical  solution  of 
the  optimal  control  problem  for  the  whole  network  system 
without  communication  constraints.  □ 

III.  Centralized  minimum  time  rendezvous 

In  this  section  we  study  the  rendezvous  problem  for  a 
robotic  network  of  first  order  agents  with  communication 
edge  map  £disk  or  Lcube  and  look  for  a  control  and  commu¬ 
nication  law  that  solves  the  problem  in  minimum  time. 

More  formally,  let  S  =  (/,  A,  Ecmm)  be  a  uniform  robotic 
network.  The  (exact)  rendezvous  task  7}ndzvs  :  Xn  — > 
{true,  false}  for  S  is  the  static  task  defined  by 

{true,  if  =  x^\ 

V (L  j)  G  Ecmm(x), 
false,  otherwise. 

for  x  =  (a;M, . . .  ,x^). 

Thus,  given  the  uniform  network  S  =  (I,  A,  Ecmm),  the 
minimum  time  rendezvous  problem  for  first  order  agents 
with  limited-range  communication  and  bounded  control 
input  is  the  following: 

minimize  X}r=o  1> 

<),P(T) 
subj.  to 

(i)  (p(-),u(-))  is  an  input-state  trajectory  of  A, 

A  =  {#},£/  =  {(Rd,U,Rd,f»ieI ,  p(0)=po; 
(ii)  i  and  j  can  communicate  if  and  only  if 
(i,j)  G  Ecmm(p^(t ),..., p[n](t))-, 

(iii)  Tmdzvs(p[1], . . .  ,pM)  =  true  for  all  t  >  T,  T  G  N. 
Here  U  is  either  f3(0,rctr)  or  C(0,rctr),  /(pW(f),  «W(f))  = 
pW(t)  -t-uM(f)  and  the  communication  edge  map  i?cmm  is 
either  L;disk  or  Ecube. 

We  refer  to  the  minimum  time  rendezvous  problem 
with  communication  edge  map  Ecmm  and  input  set  U  as 
MTlZ(Ecmm,  U). 

Next,  we  provide  some  preliminary  results  for 
the  centralized  setting  of  the  above  problem,  that 
is,  for  MTK(EcmpUU).  Let  MEB(pW  •  •  •  pN)  and 
MEO(pW  •  •  ■  p^)  the  minimal  enclosing  ball  and  orthotope 
of  points  (pM  •  •  •  pi"  1),  and  let  MBCjpM  •  •  -pini)  and 
MOC(pW  ■  •  •  pi"l)  the  centers  of  MEB(pW  •  ■  -p^)  and 
MEO(pW  ■  •  •p[nl)  respectively.  We  present  the  following 
theorem  omitting  the  proof  based  on  geometric  arguments 
because  of  space  constraints. 

Theorem  3.1:  For  all  rctr  G  R+,  p^  G  Rd,  i  G  {1, . . . ,  n} 
the  solution  of  MTlZ(Ecmpi,U),  U  =  B(0,rctT)  or  U  = 
C(0, 7"ctr),  is  not  unique  (the  problem  is  not  normal).  If 
G  B( 0,  rc tr),  i  G  {1, . . . ,  n},  then 
(i)  p(T)  =  prndzvs,  disk  =  MBC(pW(0), . . .  ,pM( 0)), 

=min{rctr,  ||pmdzvs  - 

■  vers  (prndzvs  -  (*)),  i  G  {1, . . . ,  n}, 

is  a  solution  of  M TlZ(Ecmp\,  B(0,  rctr)); 


(ii)  if  Vi  £  {1  ||pM  —  MBC(pW  •  •  •p^)||2  < 

rctr,  then  the  solution  of  MTlZ(Ecmp\,  B(0,  rctr)) 
is  given  by  p(T)  =  pmAzv s, 

disk?  Pmdzvs,  disk  ^ 
f\(E{l,...,n}  ^  ’  ^ctr)?  3.nd  —  Pmdzvs,  disk  — 

pW(f). 

Alternatively,  if  uW  £  C(0,  rctr),  i  €  {1, . . . ,  n},  then 

(iii)  p{T)  —  Prnd/vs,  cube 5  Pmdzvs,  cube  £  P-MT,  tZa^(f)  — 
min{rctr,  |pmdzvs,a  -  PaW(f)|}  Sign(pmdzvs,a  -  Pa^))’ 

i  £  {1  a  £  {l,...,d}  is  a  solution  of 

MT7Z(Ecmpi,C(0,  rctr)),  where 


P-mt  =!T 


MOC(p[1](0),...,p["](0))__(Ux_y) 


MOC(PW(0),...,pW(0))  +  -(W-ia) 


imax  is  the  largest  side  of  MEO(pW  (0), . . .  ,p^  (0)) 
and  Za  is  the  side  in  direction  a; 

(iv)  if  Vi  £  {1, . . . ,  n}  ||pM  -MOC(pW  •  •  -p^JIloo  <  Sca¬ 
the  solution  of  AfT7^.(£,cmpi,  C(0,  rctr))  is  given  by 

P(P)  —  And/ vs,  Prndzvs  G  fV(=  {  |  n }  (J  (j)  ' ,  7*ctr ) ,  and 
Ww(f)  =Pmdzvs  -pW(f)-  □ 


IV.  Distributed  consensus  on  minimal  enclosing 

BALL  AND  ORTHOTOPE 

In  the  previous  section  we  have  shown  that  minimal 
enclosing  shapes  play  a  key  role  in  the  solution  of  minimum 
time  rendezvous.  In  fact  if  the  agents  could  know  the  center 
of  such  shapes  (ball  or  orthotope)  the  solution  of  minimum 
time  rendezvous  would  be  just  a  control  law  that  drives  each 
agent  to  this  point.  Therefore,  in  this  section,  we  want  to 
explore  two  consensus  algorithms  to  compute  the  minimal 
enclosing  ball  and  the  minimal  enclosing  orthotope  of  a  set 
of  given  points  in  Rd  in  a  distributed  way. 

Here  is  an  informal  description  of  what  we  shall  refer  to 
as  the  FloodMEB  algorithm: 

[Informal  description]  Each  agent  initializes  the 
minimal  enclosing  ball  to  its  initial  position,  then, 
at  each  communication  round,  performs  the  fol¬ 
lowing  tasks:  (i)  it  acquires  from  its  neighbors 
(a  message  represented  by)  the  coordinates  of  the 
minimum  set  of  points  describing  the  boundary  of 
their  minimal  enclosing  ball  and  the  coordinates  of 
their  initial  position;  (ii)  it  computes  the  minimal 
enclosing  ball  of  the  point  set  comprised  of  its  and 
its  neighbors’  set  of  points  and  its  initial  position 
(that  it  maintains  in  memory);  (iii)  it  updates  its 
logic  variables  and  message  as  in  (i). 

Before  describing  the  algorithm  more  formally,  we  need 
to  introduce  some  notation  and  state  some  properties  of 
the  minimal  enclosing  ball.  Given  a  set  of  m  points 
{gi, . . . ,  qm}  C  Rd  in  generic  positions,  we  denote  with 
MEBbndry({gi,  ■■■,9m})  the  minimum  set  of  points  on  the 
boundary  of  MEB({qi, . . . ,  qm})  that  uniquely  identify  such 
boundary.  When  the  points  are  in  generic  position,  we  let 
MEBbndry({9i,  ■  ■  ■ ,  9m})  denote  a  minimum  set  of  points 
on  the  boundary  of  MEB({qi, . . . ,  9m})  that  identify  such 
boundary.  Moreover,  let  MBR  :  F(Rd)  — >  R  the  function 


that  associates  to  a  set  of  points  the  radius  of  the  minimal 
enclosing  ball  of  such  points. 

Lemma  4.1  (MEB  properties):  Let  (),,  a  set  of  n  points. 
The  following  statements  hold. 

(i)  there  exists  a  subset  QA  C  Qn  of  d  +  1  elements  such 
that  MEB(Qd)  =  MEB(Qn); 

(ii)  for  all  Qni,Qn2  C  Qn  with  Qni  C  Qn2,  then 
MBR(Qni)  <  MBR(Qn2); 

(iii)  if  MBR(Qm)  =  MBR(Q„2),  then  MBC(Qni)  = 
MBC(Q„2); 

(iv)  the  number  of  possible  values  of  MBR(Qni).  for  all 

Qm  C  Qn,  is  finite.  □ 

Remark  4.2:  An  important  implication  of  LemmaPl.  Hi  )  is 
that  MEBb„dry({9i,  •  •  ■ ,  9™})  has  at  most  d  +  1  points,  then 
the  number  of  packets  in  the  message  sent  and  stored  by 
each  agent  is  at  most  d  +  1  and  does  not  depend  on  n.  □ 

The  algorithm  is  described  formally  in  the  following  table. 


Name: 

FloodMEB  algorithm. 

Goal: 

Solve  the  problem  of  computing  min- 

imal  enclosing  ball  of  a  set  of  points. 

Logic  state: 

W[il  =  (P[i]bndry,Po]); 

Msg  function: 

msg(a;W ,  ,  /)  =  w  M ; 

Initialization: 

P[lIb„dry(0)  =  {pW(O)}, 

Po 1  (0)  =pw(0)- 

For  i  £  {1, . . . , 

n},  agent  i  executes  at  each  time  t  £  N: 

1:  acquire  w  ^ 

(t),  j  G  N(i) 

2:  compute 

-PWbndry(t  +  1)  =  MEBbndiy  (WV(j)  (t)), 

wv(i)(t)  = 

{wUft)  j  £  A f(i)  U  z} 

3:  update  W^(t  +  1)  =  (PWbndry(t  +  l),Po\t)) 

Remark  4.3:  For  the  algorithm  to  converge  it  is  important 
that  each  agent  keeps  in  memory  the  coordinates  of  its  initial 
position  and  thus  computes  the  minimal  enclosing  ball  on  the 
points  received  from  its  neighbors  together  with  the  point 
located  in  its  initial  position.  In  fact  a  point  on  the  boundary 
on  the  minimal  enclosing  ball  of  n  \  points  is  not  ensured  to 
be  on  the  boundary  of  the  ball  of  77,2  <  n\  points.  This  means 
that  the  coordinates  of  the  agents  on  the  boundary  could  be 
taken  out  from  the  logic  variables  during  the  first  iterations. 
This  does  not  happen,  for  example,  for  the  minimal  enclosing 
orthotope.  The  result  is  a  simplified  consensus  algorithm.  □ 
We  are  now  ready  to  prove  the  algorithm’s  correctness. 
Theorem  4.4  (Correct  MEB  computation):  Let  S  be  a 
robotic  network  such  that  the  agents  can  communicate  ac¬ 
cording  to  some  communication  edge  map  Ecmm .  For  any  CC 
such  that  the  graph  remains  connected  along  the  evolution, 
the  FloodMEB  algorithm  achieves  consensus  on  minimal 
enclosing  ball.  □ 

Proof:  In  order  to  prove  correctness  of  the  algorithm, 
observe,  first  of  all,  that  each  law  at  every  node  converges 
in  a  finite  number  of  steps.  In  fact,  using  Lemma  14711  each 


sequence  corresponds  to  a  ball  whose  radius  is  monotone 
nondecreasing,  upper  bounded  and  can  assume  a  finite  num¬ 
ber  of  values.  Then  we  proceed  by  contradiction  to  prove 
that  all  the  laws  converge  to  the  same  ball  (same  radius  and 
center)  and  that  it  is  exactly  the  minimal  enclosing  ball  of  the 
n  points.  Suppose  that  the  algorithm  converges  on  different 
balls  (different  radius  or  different  center)  for  different  agents. 
Then  there  must  exist  two  agents  that  are  neighbors  and  have 
different  logic  variables  (corresponding  to  different  balls). 
But  this  means  that,  at  the  following  time  instant,  they  have 
to  compute  the  minimal  enclosing  ball  of  a  larger  set  of 
points,  then  either  one  of  them  will  take  the  value  of  the 
other  or  both  of  them  will  change  their  value  and  take  a 
common  one.  Iterating  this  argument  we  obtain  that  all  the 
agents  must  converge  to  a  common  value.  Now,  the  ball  at 
each  node  contains,  by  construction,  the  initial  position  of 
that  node.  Since  the  ball  is  the  same  for  each  node,  it  contains 
all  the  initial  positions,  then  it  is  the  minimal  enclosing  ball 
of  the  initial  positions.  ■ 

Remark  4.5:  If  we  admit  that  the  agents  have  different 
priority,  the  initial  positions  of  the  agents  can  be  shared  by 
all  the  agents  in  time  of  order  0(n2).  The  algorithm  is  the 
following.  Each  agent  sends  the  position  of  the  agent  with 
higher  (or  equivalently  lower)  priority  that  he  has  in  memory. 
Each  position  takes  time  (-)  (n)  to  spread  in  the  network, 
therefore  the  total  time  complexity  is  0(n2).  Even  if  we  did 
not  provide  any  bound  for  the  time  complexity  of  FloodMEB 
algorithm,  however  simulations  suggest  that  it  should  be  of 
order  0(n).  Moreover,  while  the  algorithm  for  sharing  the 
initial  position  needs  to  store  a  number  of  packets  of  order 
0(n),  the  FloodMEB  algorithm  needs  to  store  only  d  +  2 
packets.  □ 

Here  is  an  informal  description  of  what  we  shall  refer  to 
as  the  FloodMEO  algorithm: 

[Informal  description]  Each  agent  initializes  the 
minimal  enclosing  orthotope  to  its  initial  position, 
then,  at  each  communication  round,  performs  the 
following  tasks:  (i)  it  acquires  from  its  neigh¬ 
bors  a  message  represented  by  the  coordinates  of 
their  current  minimal  enclosing  orthotope;  (ii)  it 
computes  the  minimal  of  its  and  its  neighbors’ 
enclosing  orthotopes;  (iii)  it  stores  as  new  message 
the  coordinates  of  the  minimal  enclosing  orthotope 
computed  at  the  previous  step. 

A  more  formal  description  of  the  algorithm  is  given  in  the 
following  table. 

In  the  following  theorem  we  prove  the  correctness  of  this 
algorithm,  together  with  the  fact  that  it  reaches  consensus  in 
minimum  time. 

Theorem  4.6  (Correct  MIX)  computation):  Let  5  be  a 
robotic  network  such  that  the  agents  can  communicate  ac¬ 
cording  to  some  communication  edge  map  Ecmm.  For  any  CC 
such  that  the  graph  remains  connected  along  the  evolution, 
then  the  FloodMEO  algorithm  achieves  consensus  on  mini¬ 
mal  enclosing  orthotope.  Moreover,  it  achieves  consensus  in 
minimum  number  of  communication  rounds  given  by 

A FloodMEO  max  max  {dist^(i,  dist^(i,  zmaxa)} 


Name: 

FloodMEO  algorithm. 

Goal: 

Solve  the  problem  of  computing  min¬ 
imal  enclosing  orthotope  of  a  set  of 
points. 

Logic  state: 

_  r  /  [i]  [i]  \  i 

Msg  function: 

msg(a;W ,  ,  i)  =  tuM 

Initialization: 

p£L,«(°) =p«](o). 

Pmax,a( 0)  =  Pa' (0) 

For  i  £  {1, . . . ,  n},  agent  i  executes  at  each  time  t  £  N: 
1:  acquire  j  £  A f[i) 

2:  compute  Va  £  {1, . . . ,  d} 

Pmiri'aQ  +  X)  =  {Pnin,a  (*)} 

p\nax,a(t  H"~  1)  =  {plnax 

3:  update  Va  £  {l,...,d} 

Wla](t  +  1)  =  ( V£L>a(t  +  1  ),Pmax,a{t  +  !)) 


where  imin,a  and  imaXja  are  the  agents  that  characterize  the 
boundary  of  the  orthotope  in  direction  a  and  minimize  the 
topological  distance  from  i. 

The  time  complexity  of  the  algorithm  is  of  order  0(n).D 
Proof:  In  order  to  prove  the  correctness  and  the  time 
complexity  of  the  algorithm  described  in  Table  ??,  we  need 
to  prove  that  it  is  equivalent  to  2d  FloodMax  algorithms  for 
leader  election  (two  for  each  direction)  running  simultane¬ 
ously.  Once  we  have  proven  that,  the  results  on  correctness 
and  time  complexity  follow  from  Chapter  4  in  [1]. 

The  algorithm  is  clearly  a  set  of  FloodMax  algorithms  for 
leader  election.  In  fact  the  boundary  of  the  orthotope  in  each 
direction  a  is  given  by  the  coordinates  of  the  points  on  such 
boundary  which  are  characterized  by  the  property  of  having 
the  maximum  and  minimum  value  of  the  ath  coordinate 
respectively. 

In  order  to  prove  that  the  exact  number  of  communica¬ 
tion  rounds  needed  is  Teioocimeo,  simply  observe  that  it  is 
exactly  the  minimum  time  for  all  the  leaders  to  propagate 
their  information  through  all  the  network.  Hence  this  is  the 
minimum  time  for  every  possible  consensus  algorithm  to 
converge.  But  this  is  exactly  the  time  taken  by  2d  FloodMax 
algorithms  running  simultaneously  and  therefore  the  time 
taken  by  FloodMEO.  ■ 

In  the  following  lemma  we  give,  for  both  FloodMEB  and 
FloodMEO  algorithms,  a  bound  on  the  time  needed  by  each 
agent  to  decide  that  the  algorithm  has  reached  consensus. 

Lemma  4.7  (Termination  condition ):  Consider  a  network 
S,  where  the  FloodMEB  ( FloodMEO )  algorithm  is  running. 
Each  agent  can  decide  that  the  algorithm  has  reached  con¬ 
sensus  if  the  value  of  its  MEB  (MEO)  has  not  changed  after 
diamg  communication  rounds.  □ 

Proof:  In  order  to  prove  the  claim  we  proceed  by  con¬ 
tradiction.  Suppose  that  after  diamg  communication  rounds 


the  MEB  (MEO)  of  agent  i  (for  some  i  £  {1, . . . ,  n })  has 
not  changed  and  the  algorithm  has  not  converged  yet.  Then 
there  will  exist  a  T  >  diamg  such  that  the  MEB  (MEO)  of 
agent  i  will  change  to  a  new  value.  But  this  means  that  the 
new  value,  stored  T  rounds  before  by  some  other  agent  j, 
took  a  number  of  communication  rounds  greater  than  diamp 
to  arrive  from  j  to  i  and  this  contradicts  the  definition  of 
diameter  of  Q.  ■ 

V.  Constant  factor  approximation  of  minimum 

TIME  RENDEZVOUS  CONTROL  AND  COMMUNICATION  LAW 


Proof:  By  the  connectivity  arguments  done  before 
and  by  Theorem  14.41  and  Theorem  14.61  we  know  that  there 
exists  T  £  N  such  that  for  t  =  T  the  network  is 
connected  and  all  the  agents  have  reached  consensus  on 

MBC(pW(0),...,pW(o))  (or  MOC(pW(0), . . .  ,pw(0))). 

Since  this  instant  all  the  agents  can  move  toward  the  same 
point  (at  maximum  speed)  without  enforcing  connectivity 
constraint  anymore.  Thus,  they  can  converge  to  the  ren¬ 
dezvous  point  which  is  exactly  MBC(pW  (0), . . .  (0)) 

(or  MOC(pW(0), . . .  ,p[nl(0))).  ■ 


The  centralize  solution  for  minimum  time  rendezvous  and 
the  consensus  algorithms  studied  in  the  previous  section 
suggest  a  dynamic  control  and  communication  law  that  plays 
a  key  role  in  the  minimum  time  rendezvous  problem. 

Here  is  an  informal  description  of  what  we  shall  refer  to  as 
the  move-toward- MBC  (MOC)  control  and  communication 
law,  CCmeb  (CCmeo): 

Each  agent  initializes  its  logic  variables  to  its 
initial  position,  then,  at  each  communication  round, 
performs  the  following  tasks:  (i)  it  acquires  from  its 
neighbors  a  message  given  by  their  logic  variables 
and  positions;  (ii)  it  runs,  as  state  transition  func¬ 
tion,  the  FloodMEB(MEO)  algorithm;  (iii)  it  moves 
toward  the  center  of  the  current  ball  (orthotope) 
while  maintaining  connectivity. 

Next,  we  formally  define  the  law  as  follows.  First  we 
assume  that  each  agent  operates  according  to  the  stan¬ 
dard  message-generation  function,  that  is  msg(a;M ,  ,  i)  = 

(a;W,tu^).  Second,  before  the  FloodMEB  (or  FloodMEO) 
algorithm  reach  consensus,  connectivity  is  maintained  by 
restricting  the  allowable  motion  of  each  agent  in  some 
appropriate  manner.  The  exact  algorithm  can  be  found  for 
example  in  [6]. 

The  state  transition  function  implements  the  FloodMEB 
and  FloodMEO  algorithms  respectively,  with  logic  variables 
as  defined  in  the  two  tables  above. 

Define  the  control  function  ctl  :  x  x  Ln  — >  for 

each  agent  i  £  {1, . . . ,  n}  by: 

ctl =  max{A*  •  {prndzvs(w[l] -pW),rctr} 

■vers  (Prndzvs{w^  -pH), 

with 

Prndzvs(wW,yW)  =  MBC(«;M ,  yW)  (2) 

and  A*  is  chosen  in  order  to  maintain  connectivity  until 
consensus  is  reached. 

In  a  network  with  communication  edge  map  Ecmm  =  Ecu be 
the  procedure  described  above  is  applied  separately  in  every 
direction  a  £  {1, . . . ,  d}. 

The  correctness  of  the  two  control  and  communication 
laws  is  proven  in  the  following  lemma. 

Lemma  5.1  (Correctness  of  CCmeb  and  CCmeo)-'  On  the 
network  S  with  communication  edge  map  T’h;ln  or  Ecu be 
and  bound  on  the  ith  control  input  uW  £  H(0,rctr)  or 
mW  £  C( 0,  rc tr),  the  control  and  communication  laws  CCmeb 
and  CCmeo  achieve  rendezvous  at  MBC(pW  (0), . . .  ,p^n  1  (0)) 
and  MOC(pW(0), . . .  ,p["l(0))  respectively.  □ 


A.  Time  complexity  of  CCmeb  and  CCmeo 

In  the  previous  lemma  we  have  proven  that  the  control  and 
communication  laws  CCmeb  and  CCmeo  achieve  consensus. 
Now  we  ask  how  fast  these  laws  are  depending  on  the  control 
bound  rctr  and  the  number  of  agents. 

Theorem  5.2:  For  rcmm  £  R+,  d  £  N,  consider  the 
network  S  with  communication  edge  map  either  Pdisk  or 
-Ecube-  The  following  statements  hold: 

(i)  for  itM  £  B(0,rctr),  i  £  {1  ,...,n},  the  control  and 
communication  law  CCmeb  asymptotically  converges 
to  the  minimum  time  rendezvous  centralized  solution 
MTlZ(Ecmp\,  B(0,  rctr))  as  rctr  — >  0+  (for  all  fixed  n). 

(ii)  for  1  £  C(0,rctr),  i  £  {l,...,n},  the  con¬ 

trol  and  communication  law  CCmeo  converges  to 
the  minimum  time  rendezvous  centralized  solution 
MTlZ(Ecmpi,C(0,rca))  for  rctr  — ►  0+  (for  all  fixed 
n).  Moreover,  it  is  a  constant  factor  approximation 
of  MTlZ(Ecmp\,  C(0,  rctr)),  he.,  TC(7(11dzvs,  CCmeo)  € 
0(^-)  for  rctr  — *  0+  and  n  — >  +00.  □ 

Before  proving  the  theorem  let  us  state  a  useful  lemma. 

Lemma  5.3  ([8]):  For  all  pointsets  Pi  C  P,  we  have 
MBC(Pi)  £  MEB(P)  and  MOC(Pi)  £  MEO(P).  □ 
Proof:  [Theorem  15.21  The  line  of  proof  of  the  two 
statements  is  the  same,  hence  we  prove  only  the  second 
one  which  is  stronger  due  to  the  stronger  result  on  the  time 
complexity  of  the  FloodMEO  algorithm.  Using  the  previous 
lemma  we  know  that  for  all  t  £  No,  then  p^hdzvs (t)  £ 
MEO(pW(0), . . .  ,p[”l(0)),  where  p^ndzvs  is  defined  as  in 
0.  This  implies  that,  once  the  consensus  is  reached,  the 
time  to  rendezvous  is  upper  bounded  by  the  time  of  the 
centralized  solution.  Hence  the  following  bound  on  the  time 
of  convergence  of  CCmeo,  Tmeo,  holds: 


Tmeo  < 


diam(pW(0), . . .  ,p^(0))  •  — 


Tctr 


TfIoocIMEO- 


(3) 

The  first  statement  is  proven  by  observing  that  TpioodMEO  does 
not  depend  on  rctr,  therefore,  as  rctr  — ►  0+,  I'weo  converges 
to  the  optimal  value  of  the  centralized  case. 

In  order  to  prove  the  second  statement,  observe  that 
diam(pW(0), . . .  ,p[”l(0))  <  (n  -  l)rcmm  and  TFloodMEO  £ 
0(n).  The  result  follows  by  substituting  these  bounds  in  0. 


Remark  5.4:  The  previous  theorem  confirms  the  intuitive 
idea  that,  if  the  communication  is  much  faster  than  the  mo¬ 
tion  (rctr  small),  then  the  optimal  solution  in  the  distributed 
case  converges  to  the  one  of  the  centralized  case.  □ 


B.  Distributed  minimum  time  rendezvous  in  one  dimension 

In  one  dimension  (all  the  agents  spread  on  a  line),  we  can 
find  a  condition  on  rctl  ensuring  that  the  move-toward-MBC 
algorithm  is  the  solution  of  B(0,  rctr)). 

Theorem  5.5:  For  d  =  1,  let  imax  and  imjn  the  agents  in 
the  network  S  with  the  maximum  and  minimum  positions.  If 
fctr  <  \rcmm  and  both  maxje^(^){p[*m»](0)  -pw(0)}  > 
2rctr  and  maxjeA^(imin){p[-J](0)  -pllminl(0)}  >  2rctr,  then  the 
control  and  communication  law  CCmeb  solves  the  task  7)n(|zvs 
in  minimum  time  which  is  exactly  the  time  of  the  centralized 

solution,  that  is,  T*  =  dlam^ — — -  .  □ 

rctr 

Proof:  Consider  the  input  sequence  of  the  centralized 
solution  for  the  agent  (and  equivalently  for  imax).  It 
is  '«[im,nl  (f)  =  rctr  for  all  t  <  T*  —  1  and  — 

1)  =  MBC(pW(0), . . .  ,pw(0))  -  1).  Since  the 

rendezvous  time  is  bounded  by  the  time  that  imax  and  *mm 
take  to  reach  MBC(pW(0), . . .  ,p[nl(0)),  we  need  to  prove 
that,  as  long  as  the  consensus  on  the  minimal  enclosing  ball 
is  not  reached,  then  zJlminl(t)  =  1(f)  =  rctr.  Due  to 

the  symmetry  of  the  problem  we  will  give  the  proof  only 
for  /min-  It  can  be  easily  shown  that  for  all  t  >  1  such 
that  pirmi  (i)  f  p[*m“l(0)  (consensus  is  not  reached),  the 
following  holds: 

+  1)  >  -  1)  + 


It  follows: 

MBCw(t  + 1)  =  +  1)  +p[iminl(0)), 

then 

MBCW(i  +  1)  >  ^  1)  +  W  +p[i-J(0)) 

=  MBC[i^(f-l)  +  ircmm 
>  MBC[im‘”](f)  -  rctr  +  ircmm. 

This  leads  to 

1  ,  1 
^"ctr  2  ^cmm  -'>  Tctr ,  ^HCl  Tc tr  ^  ^cmm  • 

The  other  two  assumptions  ensure  the  condition  for  t  =  0. 


C.  Simulations 

In  order  to  illustrate  the  performance  of  our  rendezvous 
algorithms,  we  implemented  the  move-toward-MBC  algo¬ 
rithm,  based  on  the  FloodMEB  consensus  algorithm.  We 
implemented  it  in  the  plane,  d  =  2,  over  the  disk  graph.  The 
simulation  run  is  illustrated  in  Figure^  The  32  agents  have  a 
bound  on  the  control  inputs  rctr  =  0.1,  and  a  communication 
radius  rcmm  =  3.  The  initial  positions  of  the  agents  were 
randomly  generated  over  the  rectangle  [—6,6]  x  [—3,3], 

The  FloodMEB  law  converges  in  five  steps,  while  the 
rendezvous  is  achieved  at  T  =  58.  As  it  clearly  appears 
in  the  figure,  once  the  consensus  on  the  minimal  enclosing 
ball  is  reached,  all  the  agents  move  toward  the  center. 


Fig.  1.  Evolution  of  the  network  (in  filled  blue)  according  to  CCmeb  with 
evolution  of  FloodMEB  (green  circles  connected  by  dashed  red  line) 


VI.  Conclusions 

We  have  presented  some  simple  algorithms  on  how  to 
compute  optimal  enclosing  shapes  for  pointsets  via  dis¬ 
tributed  computation.  These  algorithms  are  then  used  to  pro¬ 
vide  efficient  solutions  to  distributed  rendezvous  problems 
for  synchronous  robotic  networks.  For  future  work  we  en¬ 
vision  characterizing  the  time  complexity  of  the  FloodMEB 
algorithm  and,  in  turn,  of  the  move-toward-MBC  control  and 
communication  law. 
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